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Abstract 


In present years, the increasing use of power electronic based nonlinear devices in various power system 
applications causes harmonic distortions in source side voltages and currents. The unwanted harmonic 
distortion can introduce detrimental effects on electrical equipment and temperature increases in neutral 
conductors and distribution transformers. In this paper, a new optimization method of Meta-heuristic 
firefly algorithm is proposed for the optimal design of passive harmonic filters to mitigate the harmonics 
existing in the power system environment. The performance of the proposed firefly algorithm based 
passive filter evidences better result compared with the other algorithm based methods. 

Keywords: Harmonics, THD, Firefly Algorithm, Nonlinear loads, Point of Common Coupling 


1. Introduction 


Harmonic distortions are one of the important 
power quality issues. Due to the widespread use of 
non-linear devices, significant amount of 
harmonics present in power system. When 
distorted harmonic currents interact with the 
impedance of the circuit, voltage also highly 
distorted at the point of common coupling (PCC). 
The following problems may incurred due to 
distorted current and voltage such as overheating 
of electrical equipments, false tripping of 
protection devices and overheating of capacitors 
for power factor correction. However, harmonic 
mitigation schemes are necessary to be considered 
for maintaining good quality and reliable power 
supply. IEEE Standard 519-1992 or IEEE 519- 
2014 (1993, 2014) [1-2] recommends a guideline 
to limit and mitigate the harmonics present in 
power system. Harmonic filters such as Active 
filters, passive filters, and hybrid filters are 
recommended to reduce the harmonics within the 
limits as per the guidelines of IEEE standards. 
Active harmonic filter is one of the main filter 
techniques which consist of active components to 
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absorb harmful harmonics in power system. The 
drawbacks of active filters are more complex and 
expensive method compared with other methods 
(Subjak, Mcquilkinl990 & Bosel992) [3-4]. 
Passive harmonic filters are widely used for 
harmonic minimization studies because of the 
following advantages such as simple construction, 
easy to install, lower cost and less maintenance 
(Kazemi-Robati, Sepasian, 2019)[5]. In general, 
the passive harmonic filter provides better 
performance characteristics which divert the 
unwanted harmonic currents near its tuned 
frequency. Further, the harmonic distortions in 
voltage also reduced. In general, the conventional 
design technique for passive harmonic filter is a 
simple trial and error method but it has certain 
problems as they do not provide optimal results in 
several occasions. Nowadays, many optimization 
algorithms are developed for optimal selection of 
harmonic filter parameters.A swarm-based 
Improved Artificial Bee Colony (ABC) algorithm 
is proposed to design passive harmonic filter 
parameters (Bose,Biswas, Vasilakos,Laha 
2014)[6]. (Murugan, Rajaram and Renukadevi 
2016)[7] proposed a harmonic suppression method 
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by installation of Dragonfly Algorithm (DA) based 
passive harmonic filter in the circuit. A solution 
technology based on an optimization method of 
Genetic Algorithm (GA) is applied to determine 
the filter parameters (Maheswaran, Rajasekar, 
Ashok Kumar 2014&Yaow-Ming Chen 2003)[8- 
9]. Many other design techniques (Yang 2018) 
[10] such as simulated annealing, genetic 
algorithms and differential evolution methods are 
complex methods but have problems in 
convergence. In this work, particle swarm 
optimization algorithm and firefly algorithm are 
used to design the optimal parameters of passive 
filter. The numerical results in the table shows 
that the comparison between PSO and firefly 
algorithm and also indicates that the superiority of 
the algorithm for the design of harmonic filter. The 
remaining sections of this paper are well structured 
as follows. Section II briefs the basics of proposed 
algorithms of FA and PSO for the filter parameter 
selection. Section III indicates the performance 
characteristics of LC passive filter and 
thesimulation results.A conclusion of the work is 
presented in section IV. 


2. Proposed Algorithm 
2.1 Particle Swarm Optimization 


Particle swarm optimization (PSO) is an 
efficient and faster metaheuristic algorithm first 
recommended by Kennedy and Eberhart in 1995 
(Kennedy and Eberhart 1995) [11]. PSO algorithm 
is a self-correctedpopulation-oriented optimization 
method which is motivated by the movement of 
organisms such as bird flocking and fish schooling. 
An individual in a swarm also called a particle 
looking for optimal position in the problem search 
area and each particle fine-tunes its position from 
its present velocity, previous knowledge and the 
experience gathered by the adjacent particles. Thus 
the particle swarm algorithm can be used to solve 
complicated problems. PSO algorithm is designed 
with anassembly of random particles (solutions) 
and then examines for targets by modernizing 
generations. For each and every computation count, 
each particle is updated with two "best" fitness 
solutions. The mainand first one is the best fitness 
solution and also that fitness value is stored. This 
best fitness solution value is named as Ppeg and 
another "best" assessment is identified which also 
attained by aparticle element in the population. One 
best fitness solution is named as global best (pest) 
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and another finest fitness solution is a local best 
(best). In a PSO algorithm, the following equations 
are considered for updating of every particle’s 
significant characters such as velocity and position 
in the search space. Hence the movements of 
particle through a search space can be measured. 


Viatm+1) = w X Vigim)+c,xX ran; (piq (m)— 


Xjq(m))+C2X rano(Pea(m) —xjq(m)) (1) 
Xia(m+ 1 )=Xia(k)+ Via(k) (2) 
where, 


* m denotes the i" iteration 

¢ wis the inertia weight 

¢ Xiqrepresents particle position 

¢ Viq represents particle velocity 

* Piarepresents best position 

* (j,C2 represents cognitive and social parameters, 
¢ ran), rans are random values between numbers 0 
and | 

2.2Firefly Algorithm 

Nature inspired Firefly Algorithm (FA) is a 

meta-heuristic optimization algorithm was first 
proposed by Xin-She Yang [10] in the year 2007 
which is stimulated by the behaviour and light 
flashing arrangement of fireflies. The flashing light 
isformed by the chemical process taken inside of 
the bodies of fire flies which is_ termed 
bioluminescene and produced from the lower 
abdomen. Generally, the fireflies use the behaviour 
of flash light for signal system to invite mates or 
prey and also’ this behaviour’ used for 
protectingfrom predators. 
In general, firefly algorithm is used in 
numerousarea of optimization problems for 
determining the optimal result to make something 
promising minimizing or maximizing _ the 
applications, firefly algorithm is used to solve 
economic dispatch problems (Apostolopoulos and 
Viachos 2011) [12]. 

The formulation of flashing light is associated 
with the objective function of the problem. There 
are two important disputes are considered in firefly 
algorithm that are the changing of light 
concentration and provision of attraction. The light 
concentration of fireflies I(r) is changed according 
to the Inverse Square law which is denoted (Johari, 
Zain, Mustaffa and Udin 2013) [13] as 


I(r)= an 
4 (3) 
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For a specified medium, the intensity of light I 
is well-defined with the distance r and another term 
named as stable light absorption coefficient (7), 


_ 75-1 
l=TLe (4) 
where 


e Is represents the intensity of light at the source. 
e I, indicates the actual intensity of light. 

The attractiveness or brightness of the firefly is 
proportional to the light intensity between 
neighbouring fireflies and the changing of 
attractiveness 2 is defined with the distance 
between fireflies (r) by 


p= Ber (5) 


Where Bo denotes the brightness constant at 
distance r = Q. The distance (1;;) traced between any 
two fireflies of 1 and j are residence at x; and x}. 
These are represented as a Cartesian or Euclidean 
distance using the equation (4) 


[_d 
iE an (x, 7 iar (6) 


The following key equation is utilized to 
determine the position of the particle x; 


oS ee ii( t ') t 
x, =x, + hye X, — X; J+ OE; 


(7) 
3. Results and Discussion 


3.1 Proposed System 

Figure 1 indicates the simulation model of the 
proposed test system. The test circuit comprises of 
a three-phase 400 volt, 50 Hz Supply system with 
nonlinear load of rectifier circuit 





Fig.1 Simulink model of the proposed test 
system 


3.2 Current waveform without Filter 

The performance of the proposed test circuit 1s 
simulated as shown in Fig.l. The current and 
voltage are measured at Point of Common 
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Coupling (PCC).During the analysis, it is observed 
thatthe 5" order harmonics is dominated in all the 
three phases A, B and C and also noted that the 
current harmonic distortion (THD)) in three phases 
of A, B and C respectively are 13.34 %, 14.96 % 
and 16.40 %. 


Selected signal: 5 cycles. FFT window (in red): 1 cycles 
2 = 
0 4 


2 4 





0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
Time (s) 


Fig. 2(a) Current waveform without Filter 


It is to be noted that, the obtained THD, values 
are not satisfied the harmonic limit specifications 
of IEEE 519-1992 harmonic standard. These higher 
amounts of harmonic distortions will increase 
power quality issues in the system. 


Fundamental (50Hz) = 3.053 , THD= 13.34% 


Mag (% of Fundamental) 





Fig.2(b) FFT analysis of current waveform 
without Filter 
Figure 2(a) shows the input current waveform 
without passive filter.Figure 2(b) shows the FFT 
spectrum analysis of input current waveform 
without passive harmonic filter. 


3.3 Current waveform with filter 

The passive harmonic filter circuit is designed to 
reduce the current harmonics’ using _ the 
conventional manual calculation method. The 
measured THDyvalue is 11.23 % in phase A. Figure 
3(a) shows the input current waveform with 
passive filter by conventional method.Figure 3(b) 
shows the FFT analysis of the input current 
waveform with passive harmonic filter by 
conventional method. 


WIyial ky ary ew 


@) Dsplay selected signal Display FFT window 


Selected signal: § cycles. FFT window {in red): 1 cycles 














0 0.01 0.02 0.03 0.04 0.0 0.06 0.07 0.08 0.09 01 
Time (5) 


Fig. 3(a) Current waveform with passive filter 
by manual calculation 
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Fundamental (50Hz) = 4.729 . THD= 11.23% 
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Fig.3(b) FFT Analysis of Current Waveform 
with passive filter by manual calculation 


The passive harmonic filter circuit is designed to 
reduce the current harmonics using the PSO 
algorithm. The measured THDyvalue is 9.89 % in 
phase A. Figure 4(a) shows the input current 
waveform with passive filter by using the PSO 
algorithm. Figure 4(b) shows the FFT analysis of 
the input current waveform with passive filter by 
using the PSO algorithm. 


Selected signal: 5 cycles. FFT window (in red): 1 cycles 
2 =] 
0 4 


2 
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Time (s) 





Fig. 4(a) Current waveform with passive filter 
by PSO algorithm 


Fundamental (50Hz) = 3.251 . THD= 9.89% 


Mag (% of Fundamental) 
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Fig. 4(b) FFT Analysis of Current Waveform 
with passive filter by PSO algorithm 


The performance of the firefly algorithm based 
passive filter is verified with the help of a test 
simulation circuit. The current and voltage are 
measured at PCC. The total harmonic distortion 
value of current is 4.72 measured at PCC. It is to be 
noted that the obtained THD, values are satisfied 
with the specifications of IEEE 519-1992 standard. 
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Selected signal: 6 cycles. FFT window (in red): 1 cycles 





0.01 002 0.03 0.04 0.05 0.06 007 0.08 0.09 01 
Time (5) 


Fig.5(a) Current waveform with passive filter 
by Firefly algorithm 


Figure 5(a) indicates the input current waveform 
with passive filter by using the FFA algorithm. 
Figure 5(b) indicates the FFT spectrum analysis of 
the input current waveform with passive filter by 
using the FFA algorithm. 







Fundamental (50Hz) = 1.96 . THD= 4.72% 
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Fig.5(b) FFT Analysis of Current Waveform 
with passive filter by Firefly algorithm. 


Table 1 presents the comparison of simulated 
results of L and C values for different passive filter 
circuits for phases A, B and C. 


Table.1. Determination of L & C Values 


Manual 
Calculation PSO Firefly 
Ll 
F7 0.0740 0.4500 2.0975 
PHASE (mH) 
A Cl 
0.0015 0.060 2.42977 


L2 
0.0614 0.2875 | 2.1965 


L3 
PHASE 0.0854 0.3300 | 2.0888 
" ta | SOO | em | ee 
C3(F) 0.00013 0.0075 | 2.6376 


Table 2 shows the comparison of THD values 
for different passive filter circuits in phase A. 
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Table.2. THD values measured in Phase A 
with different methods 


Without Filter THD 13.34 % 


11.23 % 


Manual 

calculation 
With Passive 
Filter THD 
Value 


9.89 % 


PSO Algorithm 


FIREFLY 


Algorithm abe 





Conclusion 


Harmonics pollution in the power system 1s 
increases due to the increasing usage of power 
electronic controlled loads. Due to the large usage 
of nonlinear loads causes the power quality 
problems arises in the power system. Therefore, 
harmonic elimination is essential to improve the 
quality of the supply. In this work, harmonics 1s 
mitigated with the help of a passive filter, whose 
design is based on PSO and Firefly technique. It is 
observed that Firefly technique helps in efficient 
reduction of power system harmonics. Firefly is 
found to be very much efficient and superior to 
PSO method of filter design in harmonic 
mitigation. The Fire fly algorithm based passive 
filter results show that the satisfactory performance 
of the filter designed by using other methods. 
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